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This   note    is   of  a    theoretical    character.      Its  results   are   not 
meant   to   lead    to    immediate  applications;   however,   the    subject 
is   closely   conne cted  w  ith  topics    that  have    arisen  in   connec- 
tion with   AMP   Study   38    on   shock  waves. 

DISTRIBUTION   LIST  ' 

No.    of 
copie  s 


2  Office  of  Executive  Secretary,  OSRD 

10       Liaison  Office,  OsRD 

1  Att:  B.  S.  Smith,  British  Admiralty  Delegation 

3  Chief,  Bureau  of  Ordnance 

1  AttrR.  J.  Seeger 

1   "   P.  Keenan 

1   "   R.  S.  Burington 

4  Aberdeen  Proving  Ground,  Ordnance  Research  Center 

1  Att:  Col.  L.  E.  Simon 

1   "  Hans  Lewy 

1   "  0.  Veblen 

1   "  S.  Chandrasekhar 

1       G.  B.  FistiaVowsky,  Chief,  Division  8,  NDRC 
Att:  H.  Be  the 

1  Warren  Weaver 

1  J.  G.  Kirkwood 

1  A.  H.  Taub 

1  J.  von  Neumann 

1  H.  Weyl 

1  T.  C.  Fry 

5  R.  Courant 

1       E.  J.  Moulton 
1       S.  S.  Wilks 
1       M.  Rees 


TABLE    OF    CONTENTS 

nage 

Summary II 

Figs.    0,    1,    2 IV 

Figs.    3,    4,    5    V 

Figs.    6,    7 VI 

Part   J.      Thermodynamics    and  the    Shock  phenomenon 1-19 

1.  Shocks  1 

2.  Thermodynami cal  assumptions  4 

> 

3.  The  fundamental  inequality  for  the  direction  Z  Z,  -   9 

4.  Entropy  and  parametrization  of  the  Hugoniot  Line  -  12 

5.  The  limit  of  pressure  along  the  Hugoniot  Line  15 

6.  The  example  of  the  "deal  gas  16 

Part  II.   Reflection  of  a  Shock  Wave  19-29 

7.  The  algebra  of  reflection 19 

8.  Why  glancing  incidence  is  impossible  24 

Part  III.   The  Problem  of  the  Shock  Layer 30-41 

9.  Formulation  of  the  problem 30 

10.   Character  of  the  two  singular  points  Z  ,  Z,  : 

lTode  and  saddle 34 


RESTRICTED 

II 


SUMMARY 


The   following    investigation   trios    to   clear    up   the 
general  hydrodynamical   and    thermodynamical   foundations   of   the 
shook   phenomenon.      The    main   results    are    contained    in    Part   I, 
.§§1-4,    and    Part  HI,   §§9-10.      The   first   part   answers    the 
question:      What    are    the    conditions   for    the    equation   of    state 
of    a   fluid    under    which   shocks   with  their   distinctive    quali- 
tative  features  may  be    produced.      Those    conditions,    enumerated 
in  §2,    are    partly  of    differential,    partly   of   global   nature* 
Port    III   investigates    the   physical   structure    of   the    shock 
layer  whose    "infinitesimal"    width   is    of   the    order   of  magnitude 
t      rovided   heat    conductivity  and    viscosity   are    small   of    the 
same    order.      Initial   state    and   final   state    are    singular   points 
for    the   differential  equations    of   the    shock  layer,    and    it    is 
shown    that    they  are    of    such   a    nature    as    to  make   one    expect    j" 
problem    to   have    a    unique    solution. 

§§5-8   are    only    loosely    connected    with    the    main    issue 
The    author's    previous    note,    ''A   scheme    for    the    computation   of 
ock  waves    In    gases  end   fluids",    dealt    with   ideal   gases   for 
which  the   specific  heat    is   not    constant*       $6  describes   how 
the    two  basic    theorems    ( H,  )    and    ( H- )    of   that    paper    follow    from 
the    present  more   direct    and  more    general  argument.      The    prob- 
lem  of   reflection    treated    in    Part    II,  §§7    end    8,    is    not    of 
primary   importance    for    a    theoretical   understanding   of   t] 
mechanism   of    shocks,    but    is   of  t    practial   importance. 
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It    is    shown  hero    that    the    algebra    of   that    problem  is   entirely 
independent    of    the    underlying  thermodynamics    (as    is    also 
abundantly   clear   from  von   Neumann's    studies),    and    this    alge- 
bra   is   brought    into   what    I   believe    to   be    its   most    satisfactory 
mathematical   form   (by  using  6"  =   tan     <5 -,    .    tan    o  2   rather   than 
t    =   tan  o,/tan   o?   as   the    parameter    which  satisfies    a    quadratic 
equation  and    serves    to   determine    the    angles     J-,    and     S0   or" 

J-  i*j 

incidence    and  reflection   in    terms    of    the    three    states    Z-,  ,    Z   , 
ZP).      One    remarkable    general   proposition   concerning    the    im- 
possibility  of   glancing   incidence    in   the    limit    for   Infinite 
pressure    ratio    is    proved   to  hold   under    our   v;ide    thermodynamics! 
assumptions    if   one    adds    one   further   quite    natural  hypothesis 
that    is   discussed    in  £5.      But    all   this    is   of    little    practical 
value    as    long   as    the    transition  from   regular   to  Mach  reflection 
is    still  a  mystery,    and    to    solve    that  mystery   it    solids    wiser 
to   follow   von   Neumann    in    limiting    the   investigation    to    the 
simplest   cases    (although   no  explanation   is    likely   to  be    accept- 
able   unless   it    carries    over    to    the   general   situation  here 
considered  ). 

This  memorandum   is    written   without    knowledge    of 
H.    A.    Bethe's    report,    "The    Theory  of    shock  waves   for   an 
arbitrary  equation   of    statu1',   Division   8,    NDRC,    OSRD   Report 
No.    545,    to  which   I  had    no   access. 
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Pig.  0  shows  the  adiabatic 
through  ZQ  and  the  Hugoniot  lino 
in  a  p,V-  diagram;  moreover,  in 
the  form  of  a  shaded  profile,  H  as 
function  of  the  distance  t  along 
a  ray:   H  reaches  its  maximum  at  Z'  , 
passes  through  zero  at  Z,  ,  and 
is  negative  at  Z°. 
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SHOCK  V.AVES  IK  ARBITRARY  FLUIDS 
by  Hermann  Weyl 

nart  I.   Thermodynamics  and  the  Shock  phenomenon 

1.   Shocks 

At  any  plaee  and  time  a  fluid  is  in  a  definite 
thermodynamical  state  Z  and  in  a  definite  state  of  motion, 
the  latter  being  characterized  by  the  components  u,  v,  w 
of  the  velocity~u  In  a  Cartesian  system  of  coordinates 
x,  y,  z.    The  possible  thermoaynam? cal  states  form  a  two- 
dimensional  manifold  <v  ;   examples  of  quantities  having 


definite  values  in  a  definite  thermodynamical  state  are 
density  P    ,  pressure  p,  temperature  T,  volume,  potential 
energy,  heat  content  and  entropy  per  mass  unit, 
V  =  9  "  ,    E,  c    =E  +  pV,  m      .    On  ■'/   the  fundamental  thermo- 
dynamical equation 
(1)        dE  =  T  d."7  -  pdV 

holds.   It  is  once  f#r  all  assumed  that  we  are  dealing 
with  a  stationary  state  (Z,u  );  i.e.,  that  neither  Z  nor 
u  vary  in  time,  but  depend  on  x,  y,  z  only. 

In  the  limit  for  vanishing  viscosity  and  heat 
conductivity,  a  discontinuity  is  possible  such  that 
(Z,  u)  has  a  constant  value  (Z  ,  u  )  for  x  <  0  and 
another  constant  value  (Z-,,  ul  )    for  x  >  0.   The  symbols 
[0],  [1]  refer  to  the  two  regions  x  <  0,  x  >  0  which  are 
divided  by  the  shock  front  x  =  0.   Because  of  conservatioa 
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of  mass,  P  u   =  fiui  »  an<^  denoting  by  M  this  common  flow 
of  mass,  we  may  write 

o      o'     1      1 
Then  the  laws  cf  conservation  of  momentum  and  energy  give 

the  further  relations 


0              0 

Mu1   +   plf 

Mv      =   Mv, , 
o             1' 

Mw      -   Mw,  , 
o             1' 

M 


L 


+  -  (u.2  +  v„2  +  wo2)j  =  M\  i1   +  \    (ux2+  vx2+  w^)} 


0 


whereas  the  law  of  increasing  entropy  ("More  entropy 
flows  in  than  out")  requires 
Mvt^  =   M7q 
The  phenomenon  that  results  if  M  f-    0  is  called 
a  shock.    The  problem  then  spli  ts  into  two  parts,  one 
Involving  only  the  normal  velocity  component  u,  the  other 
referring  to  the  tangential  components  and  stating  that 
they  go  over  the  shock  front  unchanged: 
(2) 


(2') 
(2") 

(3) 


u 


MVo,    u1 


MV]_, 


MuQ    +   po      =  Mux    +   Pl, 

1      2    _  .            1      2 

1       +    — U           =  1,      +   rrU,  c 

0           2°  121 


V 


V,  ,     W 
1'         O 


w 


1 


If,  however,  M  =  0,  we  obtain  a  vortex  sheet  or  slip 

stream   characterized  by    the  relations 

u      =  u,    =   0,       p      =   p,     : 
o  1  '       ^o        *± 

the  fluid  in  [1]  glides  tangentially  over  that  In  [0]  and 
has  the  same  pressure  on  both  sides.   Our  object  will  be 
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the  study  of  shocks  in  an  arbitrary  fluid. 

After  introduction  of  the  values  (2)  of  the 
velocities  (2')  gives 

(4)  M2  =  Pi  -  Pq 

Vo  "  Vl  ' 
and  thereupon  (21')  yields  a  relation  between  the  two 
the rmo dynamical  states  Z   and  Z, ;  namely, 

h  '   i0  =  |m2(vo  -  vx) . (VQ  +  vx)  =  l(Pl  -  Po)    (Vq  +  vx)  oi 

El  -  Eo  =  ^1  +   Po>(Vo  -  V 

Hence  the  problem  of  shocks  is  reduced  to  a  study  of  this 

relation  between  two  states 

H  =  H(ZrZo)  =  (El  -  E0)    -  |(P]  +  P0)(Vo  -r   Vl)  =  C 
(Hugoniot  equation). 

'  e  are  only  interested  in  shocks  in  one  and  two 
dimensions.   If  in  the  latter  case  the  normal  unit-vector 
of  the  shock  front  In  the  direction  [0]  — S*  [1]  is  (■:*  (3 ) 
instead  of  (1,  0)  the  formulas  (2'),  (3)  for  the  velocities 
read  as  follows: 

(5)  c/.uo  +  /3vq  =  MVQ,   ocux  +/3Vl  =  MV1, 
(5')  -/3uQ  +  ocvQ  =   -  /^ux  +  Xv1. 

A  simple  analysis'  (which  will  presently  occupy 
us  in  much  greater  detail)  reveals  that  the  value  of  M 
resulting  from  (4)  is  of  the  order   i,q  where  q  is  the 
acoustic  velocity.   Thus  we  are  dealing  with  a  situation 
like  this:   A  certain  quantity  might  be  zero,  but  if  it 
is  not,  it  is  even  2s  1.   Clearly  quite  different  circum- 

RESTRICTED 


-4*  RESTRICTED 

stances  must  be   responsible   f»r    the    two   phenomena   of 
shock  find   slip   stream   and  M  =  0    is    in  no  way   to   be    con- 
sidered here    the    limiting    case   of  a   non-vanishing  M. 
The   problem  of   the    shock   layer  will    be    compared   to    that   of 
the    slip   layer    in   §9    by    taking   viscosity   and    heat    conduct- 
ivity   into   account  and   then   ]e  tt?ng   them   tend    to    zero.      In 
particular,    this   passage    to   the    limit  will   explain  why  a 
shock   is   not    conservative   with  respect    to   entropy    though 
it   conserves   energy.- 

2.   Thermodyncmical  assumptions 

Next  we  specify  our  assumptions  concerning  the 
thermodynamical  behavior  of  our  fluid.   They  will  be  enumer- 
ated, I-IV. 

I ,   Infinitesimal  adiabatic  increase  of  pressure  effects 
compression : 

(g)-  <  °- 

II-,   The  rate  of  compression  -  -r—   diminishes  in  this 


process  :•* 


d2v      >  0. 
dp2/  ad 


These  local  hypotheses  will  be  supplemented  by  two  assumptions 

of  "global"  character. 

III.   In  the  continuous  process  of  adiabatic  compression 

one  can  raise  pressure  arbitrarily  high.. 


•;  it  can  be  made  plausible  that  condition 
II  is  essential  for  the  formation 
of  a  shock  wave* 
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IV ,   The  state  Z  is  uniquely  specified  by  pressure  and 
specific  volume,  an d  the  point s  (p,V)  representing  the 
possible  states  Z  in  a  (p,  V)  -  diagram  form  a  convex  region.. 

It  would  be  more  natural  but  less  elementary  to 
divide  Iv  into  a  local  and  a  global  part,  the  local  postu- 
late asserting  that  in  the  neighborhood  of  a  given  state  Z 
the  variables  p,V  can  be  used  as  parameters  for  the  speci- 
fication of  states;  in  other  words  that  the  projection 
Z  ->  (p,V)  of  the  manifold  j/   of  states  Z  upon  a  (p,V)  -plane 
is  locally  one-to-cne.   The  global  part  would  assert  that 
the  projection  of  4   covers  a  convex  region  /  in  the  (p,V)- 
plane  and  that  Z   being  any  given  state  and  (p  , V  )  its 
projection,  one  never  runs  against  an  obstacle  when,  start- 
ing at  Z  ,  one  lets  Z  vary  so  that  its  projection  (  ,\l)    follows 
a  given  pa.th  in  L   beginning  at  (p  ,V  )  ("continuation").- 
All  this  means  that  K     is  a  covering  surface  over  the 


convex  X    without  signularities  and  relative  boundaries. 


But  the  convex  <y  being  simply  connected,'  this  covering 
surface  necessarily  consists  of  a  single  sheet,  or  (p,V) 
determines  Z  uniquely.    '  e  thus  fall  back  on  the  more 
elementary  formula  tion  IV.    The  assumption  of  convexity 
for  the  region  4  will  prove  quite  essential  for  our  in- 
vestigation. 

Because  of  IV,  entropy  .-.  is  a  (single-valued) 
function  of  p  and  V,  n   =  >/(p,v).   The  adiabatic  process 
of  compression  Is  thus  defined  by 
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^dp  +  ^dV  =  0 

and  since  -*-2,  ~£      cannot  vanish  simultaneously,  eq.  (29), 

V.)    7)  I 

hypothesis  I  requires  that  ^r1-,   ^rh   are  of  the  same  sign: 

>  p     '  >V  7)  p     '  3  V 

As  the  region  4  is  connected,  the  one  set  or  the  other  will 
hold  everywhere .   We  assume  the  first: 


>   Q,   %-l   >   0. 


(The  other  alternative  would  make  little  difference.) 

Infinitesimal  adiabatic  compression  is  now  described  by 

dp  =  a*,dt,   dV  =  b#.dt;   a-::-  =  M,   b-::-  =  - 1~ 

(a  posi tive  dt  corresponding  to  an  increase  of  pressure) . 

Hypothesis  II  yields 

d    /   b',N\  ^  «        n„     v.-»-da'"  -::-   db"""     .      ~ 

"XT  I  ~ -  )  >  9        or      b  -rx      -a     ■3-c      <      3, 

Qt^av/  d  t  d  t  ' 

1-e*    b    (    r^   a        +    Vrrb    /  -      a     —-a---   +  ^—  b"  /  = 

*2db*         *.#/3a*    >b*^    .    h#23  a* 
(lit)   ^     /£»  XS   _        ^     ^    ^  ^     ^,  2    <  c< 

^p2(^)       3p^v?vJp     av«^p/ 

In    the    future  we    shall  use   hypotheses    I  and  11    in   their 
analytic    form    (I')    and    (II1). 

Our  next  move    consists   in  developing   a  number 
of    consequences   from  the   assumptions   I-III.      Let    (a,b)    T    (0,0) 
be    two  given  numbers  which  determine    the   direction   of  a 
straight   line   through  Z     =    (p   ,    V    ) , 

p   =   po   +   at,      V   =  Vc   +   bt, 
the   half-line    or    ray   being    obtained   by    the    restriction    t   ^>  0 
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of  the  parameter  t.    We  follow  the  straight  line  or  the 
ray  as  long  as  it  stays  within  -Y,  end  form 

dt   3p    ^V    ^p^       Sp^V     3v^ 
According  to  ( I '  ) ,  "7  is  positive  (negative)  and  hence  ^l    in- 
creasing (decreasing)  along  the  whole  straight  line,  provided 
a  2?  0,  b  .^.0  (ago,  bgO). 

Lemma  1 .     If  ^  =0  for  a  certain  vr lue  of  t,  then 

d7  <  0  for  the  same  value, 
d~t~ 

"roof..    Substitute  the  values  of  a  and  b  derived 

from  7   =  0  into  -ppp-  and  use  the  inequality  (II1). 

Lemma  2.   If  7'  <  0  for  t  =  C,  then  l'<   0  for 
t  >  0. 

Proof.   First  assume  7  <  0  for  t  =  0.   Should 
7  change  sign  as  Z  travels  along  the  ray,  it  would  have  to 
vanish  somewhere.   Suppose  this  occurs  for  the  first  time 
at  t  =  t,  .   As  7  <  C  before  t  reaches  this  value,  7  passes 
the  zero  level  at  t,  ascending;  hence  -T-r-   ^  0  for  t  =  t,  . 
But  these  two  relations  1-0,    -rrf-  2:  0  for  t  =  t,  , 
contradict  Lemma  1. 

If  y\    -   0  for  t  =  0  we  have  ^f-  <  0  for  t  =  0; 
consequently1/  becomes  negative  immediately  after  the  point 
Z  has  started  on  its  way,  and  from  then  on,  as  we  have  seen, 
'I'  must  remain  negative. 
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1  o  turn  to  III.    Starting  at  a  given  po^'nt 
Z   =  (P0>  V  )  and  raising  p  conti  nuali  y,  we  follow  the  adi- 
abatic through  Z  ;  we  thus  obtain  a  continuous  monotonely 
descending  function  V(p) ,  and  according  to  III  we  can  make 
p  travel  over  the  entire  infinite  interval  p  ^>  p  .   -Vhile 

this  happens,  the  directional  coefficient  s  =  ^  "  P°  of  the 

V@  -  V 
straight  line  joing  (po,V  )    with  the  point  (p,V)'  on  the 

adiabatic  Tf   =  °7   increases  monotonely  from  a  certain  positive 

value  m0  to  +<*&.        m  denotes  the  adiabatic  derivative  -  -*£ 

and  m   its  value  at  the  point  Z  .    Consequently,  there  is 
exactly  one  point  of  the  upper  branch  p  >  p   of  the  adiabatic 
on  the  ray  from  (p  ,  V  )  of  a  given  direction  s,  provided 
s  lies  between  m   and-o^. 

The  analytic  proof  for  these  intuitively  evident  . 
statements  runs  as  follows.    Form,  the  (non-total)  differential 

dr  =  (VQ  -  V)dp  +  (p  -  pQ)dV, 
which  vanishes  along  any  ray  through  Z  .    '  e  want  to  show 
that  along  the  adiabatic  dr  is  positive  for  a  positive  incre- 
ment dp; 

r  =  (v0  -  v)  +  (P  -  po)g  >  e. 

By  hypothesis  II 

<6>   f  =  <P  "  P„>  S  >  °    f°r  P  >  Po> 

and  since  R  vanishes  for  p  =  p   it  must  indeed  be  -ositivc 


for  p  >  pQ 
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<7)     ^/np,4n>  dp>  0. 


P^  (P  *  Pj 


0 


o;  dp^ 


We  may  write 

(8)       dr  =  (V   -  V)2.ds  ' 

and  thus  the  monotone  behavior  of  s  is  exhibited  by  the 

equation 

^  =  r/(v  -  v)2 

dp    /   o 
It  is  clear  that  s  tends  to  m   for  p  — y»  p  ,  and  combination 

of  hypothesis  III  with  the  trivial  inequality  s  >  (p  -  p  )/v 
leads  to  s  — >  -o^3   for  p  — >&", 


3.   The  fundamental  inequality  for  the  direction  Z   Z-, 

For  the  differential  of 
H(ZZ0)  =  (E  -  EQ)  -  |(p  +  Po)  (VQ  -  V) 
we  find  by  (1) 
dH  =  dE  -  |(Vo  -V)dp  +  -i(po  +  p)dV  =  Tdif  -|  j (VQ  -  V)dp  +  (p-pQ)dv|j 

(9)  dH  =  Td-n  -  ~   dr, 
and  hence  r long  any  ray 

(10)  dH  =  Td^  . 

Let  Z  ,  Z,  be  two  distinct  states  satisfying  the 
Hugoniot  equation  H(Z-,,  ZQ)  =  0.   Making  use  of  the  convexity 
of  the  region  y      (hypothesis  IV)  we  join  the  two  representing 
points  (p0V  ),  (Px^)  in  the(p  V)  -diagram  by  a  straight 
segment.   Integrating  (10)  along  it  we  obtain 

(11)  H(Z   Z  )  =   rzZlTd^7  =  0. 

Prom   this    sinple    relation  we   are    going  to   deduce 
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Theorem    1,           For   any    two  distinct    states    Z    ,    Z, 
>z o        1 

linked   by   Hugoniot's   equation   H   =   0    the    following    inequal:;  ties 

hold. 

(12)       N*   =    (p3    -    PQ)     (VQ    -    Vx)    >   0, 

(15)       (p3    -    pQ)    +  mo(V1    -   VQ)    >  0r    (p2    -   pQ)    +  m1(V1    -   VQ)    <  0 

Before   proving   the    theorem   let  us   discuss   its 
physical    significance  *      We   have    seen    that    elimination  of 
the    velocities    from   the    conditions    for   a    shock   lead    to 
the    Hugoniot   equation    and    the    inequality   N"    ^  0.      V1  e    now 
real  ize    that   wo   may    omit    the    supplement.-. ry    relation   N"    >  0 
because    it  follows,    even   in    the    sharper   form  >t"    >  0,    from 
the    Hugoniot   equation. 

The    adiabatic   derivative 

dg   =   _v2#d£  =      2; 
df  dV 

is  the  square  of  the  "acoustic  velocity"  q.   Assume  V   >  V, . 

Ther.  the  two  relations  (13)  give 

(14)     m  <  M2  =  Pl  "  po   <  * 

Vo  "  Vl      1 
or  P   q  <  |?;|  <  P^q    ,    thus  confirming  a  statement  made  in 

£1  on  the  magnitude  of  M.    In  terms  of  the  velocities  uQ,u1 

given  by 

2      Pi  "  PQ  „2       2  _  Pl    po  „2 
c       VQ  -  V1    o       1    V0  -  V1         1 

(14)  may  be  written  as 

(15)  !uol  >   qQ,   lu1!<  q1? 

relat-vuly  to  the  shock  front  the  flow  in  [0]  ; |S_  supersonic, 

in  [1]  subsonic. 

—  RESTRICTED 


"11_  RESTRICTED 


Proceeding  to  the  nroof  of  our  theorem,  let  us 
travel  from  Z  along  that  ray  which  passes  through  Z-,  and 
hence  set  in  the  two  lemmas 

a  =  px  -  pQ,   b  =  V;L  -  VQ. 
Were  a  ^  0,  b  ^  0,  then  V  would  be  positive  and  hence  7 
monotone  increasing  along  the  segment  Z  Z,  ,  which  contra- 
dicts the  relation  (11).   Consequently  this  combination  as 
well  as  the  other  a  <  0,  b  <   0  are  ruled  out,  and  therefore 
ab  <  0,  or  (12)  must  hold. 

Next  consider 

1        3p    sv 
Were  Y   g  0  for  t  =  0,  thon  "*J  v'oulc  be  negative  and  hence 
*J    would  monotonely  decrease  wh'le  Z  travels  along  the 
straight  segment  from  Z   to  Z-,  ,  in  contradiction  to  the 
equation  (11).   Therefore 

a£T)   +  b(||)   >  0, 
^P/O    \3V/o 

and  this  is  identical  with  the  first  of  the  inequalities 
(13),  from  which  the  second  follows  by  interchanging:  ZQ 

ana  Z, . 

But  our  argument  shows  much  more,   while  Z  travels 
along  the  ray  from  Zq  passing  through  Z-j_,  °]      starts  with 
a  positive  value;  because  of  (11)  it  must  change  sign  be- 
fore Z  reaches  Z,  .    Put  ')'   remains  negative  from  the  moment 
it  vanishes  for  the  first  time.   Rise  and  fall  of  H  =  H(ZZQ) 
along  the  ray  is  coupled  with  that  of 01  by  the  relation 

RESTRICTED 


RESTRICTED 


dH   =  Tdw   .        Hence  H   first  rises  monotoncly   to   a   positive 
mr.x;mum   ana    then   decreases,    on    the    descent   pas si  ng  through   0 
for    Z   =   Z-,.       This   description    shows    thr.t   His  positive    from 
Z      to   Z-,    and   negative    after.       In  particular: 

On   the   ray  from  Z      through  Z,    the   point   Z-.    is 
the    only   one    which   satisfies    the   Hugoniot   equation   H  =   C. 

According   to    (12)    either   the    inequalities 
Px   >  pQ,    VQ   >  V^    or   pQ   >  ■p1,    V.    >  VQ   hold    for    two   distinct 
states   Z    ,    Z-> ,    satisfying   the  Hugoniot   equation  H(Z-,,    Z   )    =  f , 
i:e    shall    indicate   these   altern;  tives   by  Z,    >  Z   ,    Z,    <  Z 

respectively.      By    the    Hugoniot   contour  #y(for   a   given   Z   ) 
we   understand  the    locus   of   all  po'nts    Z-,    =    (p^V-,)    for  which 


1         v^lul 


H(Z,Zj    =   0  and      Z,    >   Z    . 

1   o  1  o 


(It   is   quite    essential    for   our    argument    to   pick    this   upper_ 
:1  >  Zo' 


branch   Z-,    >   !  „.)       Our    above   result  may    then   be    stated    thus: 


P-,    "   Pn 
On   the   Hugoniot   contour   s   =  — u   lies   between  m     andO°    t 

vo  -  vx 

s  is  a  uniformizing  parameter  inasmuch  as  the  value  of  s 
specifies  uniquely  the  point  Z-,  .    But  it  must  be  borne  in 
mind  that  so  f ;  r  we  have  not  yet  proved  that  to  every  pro- 
assigned  value  of  s  >  m   there  actually  corresponds  a  point 
on  the  Hugoniot  contour;  we  only  know  thr.t  there  cannot  be 
more  than  one . 

4._  Entropy  ana  p;  rami  tr '  za  tion  of  the  Hugoniot  l'ne 

For  a  moment  let  us  return  to  the  upper  branch  ^V 

of  the  adiabatic  through  Z  ,  on  which  p  increases  monotonely 

from  p   to  +o^  and  the  directional  parameter  s  =  ^  ~  j?T  fram  m 

o 
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too».   Moving  along  it  we  have  by  (9)  dH  =  -gdr  =  --in. dp 
.'.•rid  hence  (7)  implies 

Lemma  5.    H  =  H(ZZ  )  is  negative  along  Ot  . 

As  r.  matter  of  fact,  we  even  know  that  H  ?'  s 

falling,  R  >  0,  and  falling  with  increasing  rapidity, 
~   >  0,  eq.  (6).    The  explicit  formula  giving  H(Z°Z  ) 
for  any  po^'nt  Z  on  Ly  is 

H(z%)  =  "  \  f  PVp  -  -  \  AP°(p  -  p0Hp°-  p^p  <  o. 

7Po         ^'P0 
This  lemma  is  instrumental  in  establishing  the 

following  two  proposi t ionsr 

Theorem  2.    For  rny  two  distinct  states  Z  ,  Z, 

atlof;:'ng  the  Hugoniot  equation  one  has  "i  >  V      or  fl      <  J 

according  to  whether  Z-,  >  Z   or  Z,  <  Z  » 

Theorem  5.   The  Hugoniot  contour  is  a  simple  line 
starting  from  Z   on  which  s  and  V    are  monotone  increasing 
(s  traveling  from  m   to  +<z&,    7  from  y     to  an  unknown  destin- 
ation) . 

"roof.   Assume  the  arrangement  Z-,    >   Z   for  the  two 

1     o 

points  Z  ,  Z-,  linked  by  the  Hugoniot  equation.   Draw  the  ray 

P-,  -  p^ 
from  Z   passing  through  Z-,  :  its  direction  s  =  L  ° 

Vo  "  Vl 

lies  between  m  and  +00,  and  hence  the  ray  meets  the  upper 


br 


anch  of  the  adiabatic  line  ^in  a  point  Z°  =  (p°,  V°) 


jmma   3,    H  ~  J  7 


zo 


where,    according   to   Lemma    3,    H   -J         ToVn    is   negative, 

Jo    ' 

Hence  "V  on  Its  road  along  the  ray  must  have  passed  the  point 

Z   where  it  reaches  its  maximum,  and  running  down  the  hill 

have  passed  a  point  Z-,  for  which  H  =".   Td'y  =  0  before  coming 

"  2o 
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to    Z    .       Its   value     y.    at    the   point   Z,    of    the   Hugoniot    con- 
tour  is    therefore   higher   then   its  value     ^     at   Z      (Theorem  2) . 

Ag;in,    let    s   be   any    value   between  m      and  o&  . 
Thu    ray   going   out   from   Z      in   the   direction    s   meets    ty    in 
a    po'nt   Z      where   H   is    negative.      Hence  H  must   vanish   at   a 
certain  po?'nt   Z-,    between   Z'    and   Z°    (where   "y   is   already   on 
the   down   grf.de),    and    thus   a   point  on  the   Hugoniot   contour 

is   obtained    corresponding   to    this    preassigned  value   of 

'Pi    "   P^ 
s   =  _± 2  We   know    it   is    the    only    one.      The   question 

V   -  V  ' 
o     1 

vt   ised  at  the  end  of  q  3  is  answered  in  the  affirmative, 
;nd  •  e  may  now  speak  of  the  Hugoniot  contour  as  o  simple 
line  '-.-,  =<p(s)    with  the  unif ormizing  parameter  s.. 

Because  the  point  Z,  =  <p{z)    is  thought  of  as 
varying  along  the  Hugoniot  lino,  let  us  drop  the  index  1. 
Along  ^wc  may  use  the  differential  relation  dH  =  0  or, 
according  to  (9)  and  (8), 

Td^  -  -idr  =  Td/n  -  i(VQ  -V)2ds  =  0. 
Pc  cause 

d2  _(vo  -  v)2 

ds     2T 
thus  turns  out  to  be  positive,  fj    increases  all  the  way  along 
the  Hugoniot  line  with  s.   This  gives  a  new  proof  for  the 
inequality  "}..    >  J     by  integration  along  that  line, 


(v0  -  V)2 


'1    'o  — 2r ds      (si  >   m0) 


o 

Thus  Theorem  3  is  proved,  and  Theorem  2  even  by  two  different 

methods. 
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5.   The  Unit  of  pros sure  al ong  t ho  Hugonio t  line 

The  strr.5ght  segment  Z  Z°  of  direction  s(>  n  ) 

o  v         o' 

joining   ZQ  with   a   point    7     on  Of  converges    toward    the   vertical 
ray  V  =  VQ,    p  ^  pQ   if   s    tends    to   infinity.      7t   is   therefore 
naturr.1    to  make    the   additional  hypothesis 

V.      _If    (p    ,    V    )    lies    in  A    so   docs    the    entire 
vcrticrl    ray   V   =  V    ,    p      g  p   <  o^>. 

Physically   this   means   that  when    the    unit   quantum 
of   fluid   in  the    state    (p    ,    V    )    is    sealed  up   in   a   vessel 
of   volume   V      and    then  her  ted,    its   pressure    can   be    raised 
to    •'  rbi  trarily   high   values.       I  na:"nte::n: 

Theorem  4.  Under  the  additional  hypothesis  V 
the  function  p( s) representing  pressure  as  function  of  s 
along   the   Hugoniot    line,    tends    to    infi  nity   with   s  — >  o*7. 

p    tends    to   infinity   along  the    adiabatic   line 
(hypothesis   Til)    as   well   as    'long  the   vertical  V   =  V 

(hypothesis   v) ;    the    same   behavior   along   the  Hugoniot   line 
will    be    shown   to   be    r    consequence   of    these    two    circumstan  ces , 
We    a r gu e    as   f o 1 1 o w s . 

Let  h  be   an  arbi trarily   large   positive    constant. 
The    section  p     ^  p   ;g  p      +   h   of    the   vertical   v   =   V      lies 
in    \    and   hence    the    same    is  true    for   a    cert."  in    strip  3    of 
small   width  a   , 

VQ  -S  g     ViV0,    P0SpgP0  +  h 
to   the>    left   of   it.  may   ascertain  two   positive    constants 
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c,    C   such    that     ~-  5;   c,  ^   g   C    in   S,       Consider    the    section 

0   g   t  g  h   of  the   ray   p   =   pQ   +    t,    V   =  VQ    -  7l  t .       If    the 

positive    cons tant A-  which  determines    its   direction   is 

r 

less  than  ^—  that  section  lies  in  S  and  therefore 

h 

dt      Tp       A"^v  >  ° 

provided   A<  c/G.      Previously  we   have    characterized    the 
direction   of  such   a    ray  by    the    number   s   =   !//[_    .      Hence 
^    is    on  the   ascent  along  the  whole    segment   p     <  p  g  p      +   h 
of    a   ray   of   direction    s  provided    s   >w    =  max(h/$     ,     c/c)  . 
Such  a   ray   will  meet    the   Hugoniot    line    in   a   point  0(s)    = 
(p(s),    V(s)    )    after  ^    has  begun    to   decrease,    and  hence 
the    corresponding   value   p(s)    is  necessarily   larger   than 
p     +  h.      In  other   words*      Given  any  h  >  9  vie  may  ascer- 
tain  a   number  v>>  such    that   p(s)    >p      +   h   for   s   >w. 
The   statement   of   Theorem  4   is    the    customary  abbreviation 
for  this   fact. 

6.      The    example    of   the    ideal   gas 

An   ideal  gas   is    characterized   by   the   equation   of 
state   pV  =   RT .      IVe   normalize   units    so   as    to   make    the    con- 
stant  R   =   1,      The    eroigy    is   a   monotone    increasing   function 
of    temperature   alone,    E   =   f(T),    and 


1-fiF- 


'-6JV    +    log  V. 
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One   easily   computes 


(dp  \  p   f  »+i 

"Uv/ad-  v-rr~  ' 
(d2y')      _  p    wf 

VHvz/ad        v2#f'3 


2 

where  f,  f"  stand  for  ^J,  iL£   and 

dT'    dT2 

Wf   =   f'(f    +    1)    (2f    +    1)    -   Tf" 
Condition   I   is    satisfied,    II  requires 

Wf   >  0. 
Hypotheses    IV   and  V   hold,     the    convex   region  V  being  the 


quadrant  p  >  0,    V  >  0.        The   adiabatic   line    through  Z    , 


/i 


'  ^ffidT   +   log  |     =  0, 
o  o 


is  best  parametrized  in  terms  of  T:   while  T  increases  from 
T   to  infinity,  V  diminishes  according  to  the  law 

log  V/VQ  =  "  JT   (f'/T)dT, 
and  hence  p  =  T/v  tends  to  infinity  even  faster  than 
T  (p  ;>  T/V  );  thus  statement  III  is  verified.   (it  ought 


to  be  observed  that  V  is  not  sure  to  tend  to  0  in  this  pro- 

cess;  if  the  integral  J   (f'/T)dT  converges,  there  will,  be 

o 

a  positive  limiting  value.    Nor  i s  it  sure  that  V  tends 


to  o&  in  the  opposite  direction.   Thus  there  is  good  reason 

for  using  the  upper  branch  of  the  adiabatic  and  for  not 

going  beyond  the  statement  III.) 
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In   a   previous   note    I   found   that  for  an   ideal   gas 
the   Hugoniot  equation   is    solvable  by  a    simple   extraction   of 
square   roots   if   one  uses  T   as    the   parameter   of   the   Hugoniot 
15ne  :         change    the   rectangle    with  the   Siaes  T   ,    T    into   one    of 
equal    area    the    sloes    \7n,  'ir  of   which  have    the    'ncreased   dif- 
ference 


(T   -   TQ}    +   2(f(T)    -   f(TQ)    ) 
Then   the    point  Z   =    (p,V)    determined  by 

V/Vo   =    ^</To'      VV   =    &/T>        Po/P  =    "VT>    P/Po   =    $^o 
describes   the   Hugoniot   line.        The    corresponding  values  of  V 

and    s   are 

1   ~    %    =      fTT(f,/T)dT    +    log(  50Ao), 
o 

3  =  £_I_P°     =     £,    .     $  -  T0        =        p  T   .    a 


Consequently  our  investigation  yields  the  following  mathematical 
theorem: 

The  assumptions  f '  >  0,  Wf  >  C  imply 

<16)   f-^  *   f(T0)    '   ?F^r   <   f,(T)    : 


(17)  7  -  ?o  >0. 
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While  T  thus  appears  as  a  monotone  parameter 
fo»  the  Hugoniot  line  we  know  from  our  general  theory 
(Theorem  5)  that  J    is  of  the  same  character;  hence  -?J- 

ought  to  be  positive  throughout.   Indeed,  a  simple  calcula- 
tion gives  the  expression 

dj  =   f'(T)(^  +  &Q    -  2T)  +  (fr   -  T)   =   ,(T) 

dT  t(77^)  "      f 

the   numerator   of  which    is   positive    according   to    the 
second   of   the    inequalities    (16).      Hence    (17)    is   an   element- 
ary  consequence    of    (16), 

T 

)dT. 


7-  1^  f     ^ 


O 

However,  it  remains  a  challenge  to  the  mathematician  to 
find-  for  the  theorem  (16)  a  simple  direct  proof  that  is 
not  a  mere  disguise  of  the  "shock  thermodynamics"  of  an 
ideal  gas. 

Part  II.   Reflection  of  a  Shock  '■>ave 
7.   The  algebra  of  reflection 

We  study  a  two-dimensional  fluid  covering  the 
semi-plane  y  >  0  bounded  by  the  wall  of  the  x-axis  and  two 
oblique  shock  fronts  meeting  at  the  origin.   Notations  are 
indicated  by  Figure  2.   Let  (  ^  ±  fi ±)  >     ^  2  r?)    be  th0 
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normals  in  the  directions  [1]  — >      [0],  [0]  —>   [2]  respect- 
ively, and  (u, ,0) , (u  , v  ) , (u9, 0)  be  the  velocities.   Contin- 
uity  of  the  tangential  velocity  component  across  the  shock 
front  1,  eq.  (5'),  gives 

/|u,  =Au   -  <=U  .  v   or  u.  =  -  (  Au   -Oi^J  , 
Ay  1   /  1  o     1  o     1  a  ,  /-'l  o    1  o  ' 

and  thus  (5)  changes  into 


o  "^-iV  =  M1V1 


while   M-,    satisfies    the    equation 

Mx2    =    (Pl    -   PQ)/(Vo    -    Vx) 
Similar   relations   prevail   for  the   second   front: 

^2>   ~2\   +/2Vo   =  M2Vo>       </)V    ^(^%    "^2vo^    =  M2V 
We  may   eliminate    the  velocity  components   u    ,v      and   obtain  two 

relations  between  the    two   directions,    either  by  equating    the 

values   of  u    ,v     derived  from    (/•-,,   n ' ,  )    respectively: 

(13v)  M1/^1(V0    -   Vx)    =   M2/?2(V0    -    V2)     (=   vq), 

or,    bv    equating    the   values   of  u    ,v     derived   from   (  // ,  ,  ff0) 

and   from    (A1 -p/f'g)    respectively: 

(19    ) 

<l9v>  Vo    ^1^2    "   M2°V    =     "^|(V1M1o42   -   V2M2^). 
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(The    first   equation    of    this   equivalent    second   set    is    omitted 
because   it  will   not    be  needed.) 

We    discuss    the    signs   by  means   of   the    first  elimin- 
ation.     P'or   geometric   reasons     c* cxg,    c*'       /£      -  oc     AJL, 
are   positive;    hence    by    (18    )    M,  ,    Mg   are    of   equal    si  gn,    and 
since    the   regions    [1]    and    [2]    get   interchanged   and   M-,,    Mg 
replaced  by   -M?,    -M,    if  the   direction   of    the   x-axis    is   re- 
versed,   we   may   assume   NL    >  0,    Mg   >  0.      The    lav;   of    increasing 
entropy  requires 

V    %   -    V    ^=°>      M2(?2    "    V    ^°> 

which  gives  ^,  <  ^  <)  ,  and  thus  by  Theorem  2  the  arrange- 
1  —   o  =•  / 

ment  Z-,  <  ZQ  <  Zg  must  prevail, 

(20)  Pi<P0<P2'      Vl>Vo>V2 

Then    (18    )    shows  /J,    and   /9g    to  be   of   opposite    signs,    and 

thus      /S      >  0,      /^g    <  0    (because    fj ^   <  0,    /$ 2   >  0  wou  Id 

contradict    oCg  /On    -  OC   y-/g  >  0    ). 


We    "ntroduce   the    abbreviations 


j.s- 


2 


PQ   -   P]_    =   V{  ,      P2    "  P0   =   P2    ;      Vl    ~   Vo   =  Vl    *      Vo    "   V 
so    that  l\1      -   P1'"'A1'/$,      Mg      =   ?2  /V2"' 

Assuming   three    states   Z-,,    Z    ,    Zg    to  be   given    in 
agreement  with    the    inequalities    (20)    let  us    try   to    solve 
in   the    most    convenient    fashion   the    algebraic   equations    that 

RESTRICTED 


-22- 


RESTRICTED 


determine  the  normals  (^  /j^ ,     (o/g  A)  of  the  two  shock 
fronts.   The  double  expression  for  ^Q/v0   derived  from 
(18u),  (18v)  gives 

^22yo +  Av  °s  A(vo  -  v- 

This   is    a    bilinear  homogeneous   relation  with  respect   to   the 
two   sets   of   variables   t-,     :    tg  =  C*^  fl 2    :  c*2   ft^   and 

$1    :<f2   ^l^     :    A  A' 

(22)  A11t1tf"1   +   A12tl(f2   +     A21t^a    +   Ajjgtgtfg   =   0, 

with   the    coefficients 

All   =  VoV>      A12   =  V2V1*'    A21   =  VoVl*'    A22   =  VlV' 
Write    (19    )    or     CV^, (11,0c  g   -   M2°*l^    =62^V1M1^2   "  ^^^l^ 
in    the   form 

^lV^l    "   V2^2}    =   U2Ml(Vo^l    "   V1^2} 
and  multiply  at   left  and  right   by  Mg(VQ   -   Vg)A  g   = 

M-,  (V      -   V -,)/$,    respectively.      The   result    is   another   bi- 
linear  relation 

(23)  Bnt1  ri   +   B12tlS^£   +   Bgltg  ^   +   B22t2<T2   =   0 
with   the    coefficients 

Bll   =  VoP2'"'    B12   =   -V2p2*'    B21   =  Vl*'    B22   =   -VlPl'"" 
One   may   look  upon    (22)    and    (23)    as    two   linear 

homogeneous   equations    for  3~i  >6""o   by    the   elimination   of   which 

one   gets   a   quadratic   equation  for  t,  :tgj    after  solving    it 

either    (22)    or    (23)    gives  (J-[  :6~2      as   a   broken   linear   function 

of   t',  ,to.      One   may   also   work  the   scheme,    as  we    shall  do   here, 

J.         & 
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the   other   way  around,    by   first   eliminating    t.  ft~f    thus 
obtaining    the    formulas 

All<ri    +   A12«T2,       A21^    +   A22<T2 

Bll  *"l    +    B12  C2>      B21  ^1    +   B22  &~2 


=  o, 


"  Ii  =   AH   ^1    +   A12  &2    -  Bll  ^1    +   B12  ^2 


A21  €"l   +   A22 


^l^l"  +   B22r2 


(24) 


The   quadratic  equation  for  <T  :   <C,, 
l11(vq  ^)2  +  l12vo  6-x  <T2  +  l226"22  =  0 


is   readily   sees   to   have    the    coefficients 


lll=Vl  >2   "V2  Pi . »    L22   =   -V1V2(VP2*-V1  >1    ^ 

L12  =  (Pl'::"   +   P2*)(V1V2*   -   V 


*). 


We    introduce   the   non-homogeneous  variables  5*=  -   v    ey-i/Sp* 
t  =   -t,/tg,    for  which£*?is   a   permissible   value,    and  thus 
arrive    at   the   final   equations 


(21)    shows   thatCX     n     and  CX       IX  must  have 
opposite    signs,  O^   A  .  O^  /^    =  O^O^g./J,  /CL   ^  0  or  6~  2?  0, 
Vice   versa,    if    this    condition   is   satisfied,    the    second   ex- 
pression  for    t  proves   this   quantity    to  be   positive,    and 
hence   the    square   roots   of 


V 


-<&■«■'•& 


=    art 


-l 


may  be  extracted: 
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^lZ/^l   =    tan    ^1   =   /((T^o^    *°^k   =    tan  ^2   =  fa/Wo)' 

Let    the    two   states    Z      <  Z      be    given;    the   possible 

states   2g   >  ZQ   related    to   ZQ   by    the   Hugoniot   equation  may  be 
specified   by    the  directional   parameter    s   =  Pp  /Vp    .      The 
relation    (25)    then   becomes    an  equation  F(s,<T)    =   0  which  de- 
fines  a   contour    /in   a    (s,6")    -diagram.       Since    all    quantities, 
including   the    normals  of   the    shock   fronts,    are   uniquely  ex- 
pressed  in   terms   of    the  variables   s   and  6*  the    situation   is 
best  discussed  by   means    of   this    contour. 


8.      Why   glancing   Incidence    Is    impossible 

We   now   reintroduce    the    thermodynami cal   assumptions 

I-IV    of  $2      and   the    fact   that   Z,    and   Z^   as   well    as    Z      and   Z0 
>j     '  1  o  o  2 

are   linked   by  the    Hugoniot   equation.         *Ve    then   learn  from 

Theorem   1    or   eq.    (14)    that 

(27)      mx   <  P-jVv-,*   <  mQ   <  P2*/Vp*   <  m£, 

consequently  P1'"/V-L"    <  P2"/v2"    or   Lll   =  Vl*'*p2'"    "   V2'"pl'S'>  °' 


Thus  G 


1,02 


1    :    0   or  6=  C°cannot  be    a    solution   of    (24), 


2 
hence   glancing    incidence   is    impossible.      V  _^(p  -,  '  /V-,  " ) , 


2        -»-        -::- ) 
V       (P2"/V2         ar°    th°    scluaros    I  u( 


(1) 


u 


(2) 


of    the    com- 


ponents  of  the  velocity  u      in   the   direction   of    the   normals 
1   and    2.      Thus    the    inner   part   of  the    inequalities    (27)    may   be 
written  as  ,  , 


%(1)! 


<    q, 


< 


u 


(2) 


and    the  prevention   of  glancing   incidence    is    seen   to   be   due    to 

(1) 


the  separation  of  the  two  normal  components 


u 


I   (2) 
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of   the   velocity   In   the   middle  region   [0]    by   the    acoustic 
velocity   q      prevailing    there?      the    stationary    flow    In    [0] 
is    subsonic   as   far    as    the    comnonent     u         ' 

-1  I  f\ 


3  n   the    direc- 
tion  of  the  first  normal   is    concerned;    it    is    supersonic 

(2) 

with  respect   to      u        ; 

o 

We  could  have  demonstrated  the  impossibility 
of  glancing  -'ncidence  more  directly  as  follows.   The  equa- 
tion (18,  )  shows  that**,,  °*p  may  vanish  only  together,  or 
thatc*  C*  =  0  implies  the  vanishing  of  both  ^,  and  "* 

A    2  A2 

(head-on    impact).       In    that    case    /C5        =   JJ        -   1,    and    the 

second  equation    (18v)    yields  M-L2(V]    -  VQ)2   =  M22(VQ   -   Vg)2 
or  p,  "V-,  "    =  Pp"Vp".      Likewise    the    relation    (18    )    proves 
that    /J-,     Pp=©    implies  the  vanishing   of   both   /^-,    and    /<<-, 
(glancing   incidence)    and  in  that   case   we   must  have 

^l2   =C*22   =   X   and    thuS   Ml2   =   M22    0r   P1*/V1*   =  P2';:'/V2* 
by    (18    )';      It   is   the    latter   relation  which  Is  made    im- 
possible  by   our    fundamental    inequality    (27).       (By   the 
way,    these    considerations  reveal    that    the   equation    (24) 
holds   even    in   the   two   limiting    cases    of  head-on    impact   and 
glancing   incidence,    although   the   conditions  ^/g   =  ^    i\ 
-  0   then    prevailing    invalidate    its   derivation. 

I   entered    into    a  discussion   of    the   problem   of 
reflect' on   mainly   for    the   purpose    of    showing    that    the 
inequalities    (15)    for   a    linear    shock   are    responsible    for 
the   simple    fact  which    I    once    more   repeat   as 
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Theorem  5.      Un_der_  the    hypotheses    enumerated   at   T-jv 
_in  §2,    glancing   incidence    is    impossible. 

For  the    discussion   of    the    contour  J~  in   the    (s,<y) 
-diagram   it   is    convenient   first  to   ignore   the    condition^  0. 
The    contour   is    such    that    to    every   abscissa    s   >  m      there    corres- 
pond  no   or  two   values   of  the    ordinate  (5".        As   long  as   we   keep 
away  from   the    ends    s  -^  mQ   and    s  -><*=  (Twill    keep   within   finite 
bounds.      Hence,    generally   speaking,    the   contour  will    consist  of 
a   number    of   separate    lines,"    all    closed  with    the    possible   ex- 
ception of   one   opening    toward    s  — $»m  ahdone   opening   toward 
s  — ■>  &o  .      Fig.    3    shows    the   actual    situation   for   an   ideal   gas, 
Fig.    4   a   more    complicated   possibility    such    as   may    occur   for   a 
fluid    of    fairly  wild   thermodynanical  behavior.      The   horizontal 
6*  =   0    cuts   away   the    "submerged"   part    6"  <  0  which   is  without 
real    significance,    and  in    the    figure    is    indicated    by  a   broken 
line.      The    picture   makes   it   quite    clear    that   head-on    impact 
is  a   perfectly   normal    and   regular   occurrence.      Whereas    the 
details    elude   general    analysis,    we    can  deduce   from    (25)    fairly 
accurately   the   behavior   of  the    contour    toward    the    two   ends 

s  — r>  m        and    s  —=>  do. 
o 

Pi  >    V-i  >    Pn>    vn   are    fi*ed,    while   p?,    Vp   vary  on  the 
contour;    but    it  must  be   borne    in  mind    that   Vp   never  exceeds  V    . 
With   s   decreasing    to  m      the    quantities   p0",    Vo      converge 
to    zero,    and   thus    the   equation    (25)    reveals    that   from      a 
certain  moment    on  there    are    two    real   roots  G"    tending 


-::-   I   use    the    word    line    for    the    trace    of   a    continuous  motion 
(=    continuous   image    of  an    interval),     contour  for    the    locus 
of  points    in  a   plane   satisfying   an  equation    (e.g.    the 
isohypses   on   a    contour   map). 
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to  -v,  and  -  £X> respectively.   Hence  for  s  sufficiently 
close  to  its  lower  limit  m  ,  the  contour/"  has  two  branches, 
but  both  are  in  the  submerged  region  and  therefore  without 
real  significance. 

Much  more  interesting  is  the  other  limit  s  — >  00. 
Here  the  hypothesis  V  (|5)  and  its  consequence,  Theorem  4, 

Pg'  — >  ^for  s  -»  00, 
come  into  play.   Divide  the  equation  (25)  by  pp*"".   This 
done,  one  sees  that  for  sufficiently  high  s  its  coefficients 
deviate  as  little  as  one  wishes  from 

1^1  =  V*,      L^g  =  -V^gVg*,   L^g  =  V^g*  -  VgV^. 

The   discriminant   of    the  T7~   equals 

g 

~ll      -    %1T22    =    (VlV'+   V2Vl'^2    =  V^]    "   V      * 

V^l    "    V2- 

Hence  as  s  approaches  infinity  we  shall  have,  from  a  certain 

1       ff 

moment   on,    two    distinct   real    roots  6"    ,0       or   branches   of    the 

J  '• 

contour.      Lot   fi     refer   to   the   upper,  (T      to    the    lower   branch, 

r  h 

S"      >  6"    ,    and   use    the   notation   f(s)'^y'g(s)    as   an   abbreviation 

for    lim    (f(s)    -    g(s)    )    =0.        re    then    find   the    following 

S->oo 
limit   relations 

5'~v1.^l_-J/2    <r"  ~  _v       (5'  - <s"~  ^.(VUL>1 

Vl    -  Vo  V]     -   VQ 


Because  V,     .   _2 ~      is    posi  tive,       -V0    negative,    one    is 

1    V,  -  V  2 

1     o 

tempted  to  believe  that  the  upper  branch  will  finally  move 
in  the  "real"  region  <5  >   0,  the  lower  in  the  submerged  region 
(5"  <  0.   However,  such  assertions  would  be  unwarranted  sinou 
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V0   may   at    times    come    too    close    to   V      and/or   C.      Yet   a 


o 


s 


v  (v,   -  v9) 

— "     ^ — c  always    remains    greater   than  V      it   is    true    that 

1  o 

the    two  branches   never    come    close   to  one   another.      Analysis 

uses   the  formulas 

11m  f(s)    %  A,      lim  f(s)    ^  B      (for    s  -^»oo) 

to  expvess   the  following  behavior  of  a  function  f(s)  :   '  ith 

any  positive  €    one  can  associate  an  w  such  that  f(s)  <  A  +6 

or  f(s)  >  B  -£■   for  all  values  s  >  <*;.    Thus  our  relations 

imply 

: i    v  v         ' 

i  lim  6"  <,        o      lim  o  >  0  ; 

v-i-v  '  — 

(28)yLim  6"   g  0,  lim  6     ^  -  V  ; 
lim  (  £>'-  <S  )    ^  VQ. 

In  particular: 

Supplement  to  Theorem  5.   Owing  to  the  hypothesis  V, 
glancing  incidence  is  imposs i bl e  even  in  the  limit  for 
s  — B»-  c*? ,  in  the  sense  that  6" stays  bounded  while  s  approaches 
infinity. 

In  formula  (26)  the  denominator  Pp'"'(Vp  +  &  )  in  the 
"real"  region  (?**£  0  is  greater  than  Po"vo  ^  Ppvp*   ^or  idcal 
gases  it  is  known  that  PpV?  =  T?  tends  to  infinity  with 
s  — =s»  ^^  on  the  Hugoniot  line  and  that  the  temperature  T? 
may  be  used  as  parameter  on  that  line  -'n  stead  of  s.   Consider- 
ing that  CT  stays  bounded  as  Tg  ->  oC'  .  we  thus  find  t  — >   0, 
<5"t'-^»  0,  (C  — >  0,  and  therefore  the  process  of  reflection 
is  sure  to  approach  head-on  impact  (and  not  its  opposite, 
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glancing  incidence)  as  TQ  -^.  o^,    Moreover  it  is  easily 
seen  that  Ng  =  Po'Vg"  is  a  monotone  increasing  function 
of  Tp ?  and,  under  the  assumption  that  the  energy  function 
E  =  f(T)  of  the  ideal  gas  tends  to  infinity  with  T  -*»  Cx=>, 
Ng"  increases  to  infinity  with  Tp  -^>  CO ,      Hence  (25") 
reveals  that  there  is  only  one  branch  of  the  (s,S)  contour 
J   in  the  "real"  region  (ET;?  0  as  soon  as  Tp  exceeds  a  cer- 
tain limit,  namely,  as  soon  as  Np"  has  become  greater  than 
N,  .   Ml  this  refers  to  ideal  gases  only.   It  does  not 
seem  safe  to  venture  similar  statements  for  arbitrary  fluids; 
what  we  can  predict  here  is  the  behavior  of  6"  as  exhibited 
by  the  supplement  to  Theorem  5  and  the  more  precise  formulas 
(28). 

In  physical  reality  the  process  of  simple  reflec- 
tion as  discussed  in  8^7  and  8  stops  and  gives  way  to  a 
more  complicated  phenomenon  as  soon  as  s  exceeds  a  certain 
limit.   The  problems  relating  to  this  Mach  phenomenon  are 
not  yet  solved,  not  even  for  an  ideal  gas  with  constant 
specific  heat. 
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Pert   III.      The   Problem   of   the   Shock  Layer 
9.      Formula  ti  on   of   the    problem 

The   fact  that    (1)    is   a   total  differential   is 
expressed   by    the   equation 

or 

f?ql  ^T        dV       2T       ^  7  _   , 

1     j  bp  '.§v     av  •  Tp  " 

For    the   sake   of  brevity  we   denote   partial   derivatives   by 


subscripts   as   in 


2>V 


3P        'p'  av         '"" 

Our   further   argument  will    be   based    on   the    thermo- 
dynamical   assumptions    I-IV    (§2),    but   in   its    course    it  will 
become    clear    that    one   more    condition   of  highly  plausible 
nature    is    to    be    added    to   our  lists 

^-a •      Heating   a   quantum   of_  fluid   at    constant   volume 
raises    its   pressure   and    to-Tpcrature . 

It  must   be   remembered    that    the   analytic   formula- 
tion   (I')>    (II1)    °f   the   hypotheses    I  and    II  depended   on  the 
fixation   of   a    sign.      For    the   previous    investigation    this 
did  not  matter   greatly,    but   it    is    essential   for   the   behavior 
of   the    shock   layer.      We    therefore    require    explicitly   by    the 

statement    concerning  pressure    in  la,    that(-f=-rS  >  0 

V     //V   =    const. 
or     ^      >  0.        The   other  part   of   la   then   adds  the    inequality 

Tp/7p   >  0   or   Tp   >  o.         The    requirement  V   that  by   heating 
at   constant   volume    one    can  drive   pressure  up    to  arbitrarily 
high  values  will    not  be   needed, 
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In  a  stationary  field  the  state  (Z,  u)  of  a  fluid 
of  given  thermodynamic  a  1  n;.ture  depends  on  the  spatial  co- 
ordinates x,  y,  z  only,  and  not  on  time.   The  five  equations 
which,  in  the  absence  of  an  external  force,  express  conserva- 
tion of  mass,  momentum  and  energy,  take  heat  conductivity 
and  viscosity  into  account  by  the  flow  of  heat 

-fi>   =  -A  .grad  T        and  a  stress  tensor  S  of  the  form 


^ui  .^» 


+ 


k' 


S.,  =  p  d  .,  -  S,.  ,   S.i  =  ;u   divu.d..  +u.  — 
ik   ^   lk    ik'    lk  r  ik     \c>  xv   e>  x./ 

(in  the    writing  down  of  which  we   have   for  a   moment   adopted 
the   subscript  notation).     A.  -  \[Z] ,  p.,  yu '    are   given  func- 
tions   of    the  thermodynamical    state    Z ;  X  >  0,yu   >   0, 

2 

/i '    +  ->^u   >  0,  Ihey   are    small .      In    two   ways 

the   passage    to    the   limit    A  ,  yu,  /l1    — t   0    can   ;Tive    rise    to   a 
discontinuity    in   the   plane    x  =  0: 

(A)  Write  £x  for  x  and   £A,     €/\x,  €/a '    for  A,  ya,/a', 
leaving   all   other    quantities   untouched,    and    then   let   the 
positive    constant    factor  C    tend   to    zero    ( shock    layer)  . 

(B)  Let   €    tend    to   zero    after  writing    £  X,    £u    for 
x,    u    and  €2A,    C2u,    eju'    for    A,/j,  a>  '     (slip    layer)  . 

The   different  effect    is    clearly   exemplified   by 
the    continuity   equation 

Substitution    (A)    changes    it    into 
which   for  £   — *■  0  reduces   to 

|_(fu,  =  ._ 

RESTRICTED 


-32-  RESTRICTED 

whereas  the  process  (B)  leaves  the  three- term  equation  (30) 

unchanged.    This  is  typical.   None  but  the  derivatives  J*— 

survive  the  limiting  process  (A),  the  others  —-,    —   disappear 

ay      az 

from  the  equations.   The  shock  layer  is  thus  seen  to  be  s 
one-dimensional  problem,  no  interaction  taking  place  between 
the  several  one-dimensional  fibers  of  space  y  =  const., 
z  =  const.   The  slip  layer  on  the  other  hand  constitutes  a 
far  more  complicated  three-dimensional  problem,  identical 
with  that  of  prandtl's  boundary  layer.   We  are  here  concerned 
with  the  shock  layer  only. 

v.iththc  abbreviation  ac"  ~/>V    +  2/u(>   0)  the  limit- 
ing process  (A)  loads  to  the  following  equations: 
(51n)  fu    =  const.  *F  M, 


(31^)  Mu  +  p  -/"■"$$  ~  const,, 

(31  )      Mv  -/^^j  =  const.,  Mw  -/t<g£  =  const., 

(31d)   (P  -M   a;  +f  E  )u  -^(^  +  w  ^;  -/i^  -  const., 

"  12     2     2 

where  E   denotes  the  total  energy  &  +   -^(u   +  v   +  w  ). 

Here  the  quantities  A,  ^,  /U-x    are  no  longer  small.   we  assume 

that  (Z,  u)  approaches  definite  constant  values  (Z  ,  u  ), 

(Z-,,  u,  ;  for  x  — ^  -  oo  and  x  — >  +o^  respectively,  whereas 

the  derivatives  -3—  tend  to  zero. 

dx 

The  energy  equation  is  equivalent  to  the  entropy 
equation     ,  j>    Hrp 

(32) 


A1- 


W   V^W   +  \,dx7  j  +  T  (JR  )) 
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from  which   one    lecrns    that   the    flow   of    entropy 

r      '  T      dx        •    /  T  dx 

is   a  monotone    increasing  function  of    x  and  therefore   • 

M  7  ^      S>  M  *]    .      The    equality    sign   could   prevail   only    if 
u,    v,    w,    T  and  hence   by    (31,  )    also   p  were    constant.      Exclud- 
ing  this  trivial    case    of    a    constant    state      (  Z,    u)    we    thus 
obtain   the    strict   inequality 

(33)  M7  1   >  M  lo 

implying  M  =£  0.   The  latter  remark  shows  how  misleading  5t 
.is  to  link  the  alternative  M  =}=  0,  M  =  0  to  the  discrimina- 
tion of  shock  and  slip  layer;  rather  do  these  two  possibili- 
ties correspond  to  the  two  ways  (A)  and  (p)  of  passing  to 
the  limit  of  vanishing  viscosity  and  heat  conductivity. 
Being  sure  now  that  M  does  not  vanish,  we  fall  back  upon  our 
old  relations  connecting  (20,~ uQ)  and  (Z1,— u^),  including 
the  Hugoniot  relation.   Let  us  assign  the  labels  0,  1  to  the 
two  states  Z  ,  Z-,  and  accordingly  orient  the  x-axis,  so  that 
Z   <  Z-,  .   Then  Theorem  2  combined  with  the  law  of  increasing 
entropy  (33)  giv^s  M  >  0.   The  meaning  of  this  inequality  is 
that  the  shock  front  moves  in  the  direction  [1]  — >  [0],  or 


•K-Indeed,  in  general,  p  =  const.,  T  =  const,  imply  V  =  const. 
There  is,  however,  one  extremely  exceptional  case  in  which 
the  conclusion  does  not  stand,  namely  if  M  =  0  and  if  for 
some  constant  value  of  p  the  temperature  T  as  function  of  V 
stays  constant  in  an  entire  V  -interval.   Maybe  this  is  a 
warning  against  thu  danger  of  a  singularity  in  the  shock 
layer  problem  of  which  one  steers  clear  with  certainty  only 
if  one  assumes  Tv  £  0  or,  what  is  the  same,  if  the  hypotheses 
I-IV,  la  are  preceded  by  one  stating  that  p,  T  may  serve 
as  local  parameters  on  7, everywhere. 
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th    t    ZQ    and    not    Z  Initial   state    of   the   fluid    not 

yet   reached   by  the    propagating    shock. 

10 »      Character  of    the    two   singular   points    Z    ,    Z-.S      Hode 
:nd    scddle-x 

It    is   easy    to   dispose    of    the   equations    (31    )    for 
v   and   Yv.      Single    out    the    one   for   v, 

)i— -   =  M(v-b),        b=   const., 

and    introduce    the    variable 

£-  M  f*  dx       f,t(x)    =  /v[Z(x)]"l. 
(f~  rVo    ;»[x)  \  J 

■..e    find 

)    -1     =  C.ec    ,  C    =  const. 

u(::)    tends    to  ive    values    uQ    =//u[ZQ],  yu-^    =  yu,[Z-L]    for 

x  — >  -  op j      x  — >   +ofirespectively.    therefore 

£  =i;>:   +  <^(::)  for  x  ->  -°°» 
l       ™  x    +  <t(::)   for  x  ->  +C<'. 

Considering    that    M   >  0  one    realizes    that   v   cannot    tend    to 
a   finite    limit   v-^   for  x  — >   +c&  unless    C   =  0.      Thus    v  and 
us  t   be    constant    throughout    the    shock   layer. 

(31    )    is    satisfied   by  u  =  MV.      Denoting  by  a/M 
the    constant    at    the    right    of    (31,  )    and    substituting    from 
(31b)    the    value    of   j    -  /*?'  £j|  into    (31d)   we    are    left,    as    is 
readily    seen,   with   the   following    two    ordinary  differential 
equations    of   first    order  for  the    state    Z   =  Z(x): 
^<*M.g  =  h2(V   -  a)    +   •  , 


<54W  X  d: 


E.g=    :    -^2(V-a)2   -c. 


*  For   ideal   gases    of   constant  specific   heat   with  certain 
special  values    of   the    constant  y    ,    cf.    R.    Bffckor,    Ztschr. 
f.    Phys.    8    (1922),    321-362;    in   particular,    pp.    339-347. 
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The    constants    M   ,    a,    c   are    determined   by 

•'2vo    +   Po    =    :i2vi   +  Pi   =   *?*> 

Eo    "  lI/l2(Vo    "   a)2   =  El   "  !m2<V1    "   a^   =    c« 

As    the  right  members    of    (54)    vanish  for   Z  =   Z 

and   Z1 ,    the    two  points    Z   ,    Z-,    turn   out   to   be   s  ingular  points 
for   the    differential   system  (54).      If  we    care   for    the    tra- 
jectory  in   the    (p,V)    -diagram  only  along   which  Z(x)    moves 
while  x   runs    from   -  o<=>    to    +oowe    are   faced  with  one    ordinary 
differential  equation   of   first   order    in  the    (p,V)    -plane, 

.   A»cLT ^__      =  /*~  •  dV 

.    "  |m2(V    -   a)2    -    c  (V    -   a)    +    (p/M2) 

But    contrary   to   the    customary  viewpoint    taken    in   the    classical 
investigations    on    singular    points,    which  were    started  by 
Briot   and  Bouquet   almost   one   hundred  years   ago    and  are   re- 
produced  in  all   treatises   on   differential   equations,    our 
interest    is    in   the    parametrized   trajectory   Z   =   Z(x). 

Let   us   briefly  recount   the  well-known  facts   about- 
singular   points    relevant   for   our    purposes.       Placing    the 
singular   point   at    the   origin  0   we   study   two   simultaneous 
differential   equations   for    the    ;!vector;i   <p-    (^(x),     tfSU(x)    ) 
in  a    <n ,  ci£   -plane, 

-TET-lllfvW'     -E5  =  V  ?i>  ?J 

where  F,,  F?  are  given  functions  of  ^  ,  cP^   vanishing  at 

the  origin.   From  the  outset  it  is  evident  that  any  trans- 
lation x  — >  x  +  h  carries  a  solution  <P  =(  ^(x),^p(x)  ) 
into  a  solution  (''equivalent  solutions").   Suppose  the 
Taylor  expansions  of  F-, ,  F?  begin  with  the  terms 

RESTRICTED 


-  "-Y 


-36-  RESTRICTED 

pi{/rf2^Riifi +  R]2^2'  vft'jy  ~ 

R2lfl    +  R22^2* 
The   equations 

f  1   =  R11?1   +   R12?2'    ^2   =  R21?1   +   R22/2 
define    a   linear  mapping  ^->  f    =  R^   •      We   determine    its 
two   eigen-values   k  and   corresponding   ei gen-vectors   &  =  ■> 
(or   rather  ejgen-d  rections)    by  Roc  =  k.x  ,    explicitly 


(k    -    R11)<x1 


"    R12*2   =   °> 


-    R21^1   +    (k    -    R22W2    =    0, 

and  assume    that   the    two   roots  k   ,    k„   of   the    secular   equation 


k  -   R 


11' 


-   R 


12 


-    Rot  >        k    -    R0, 


=    0 


v21'       '         ll22 
are   real,    distinct  and   ^  0.      In  a    suitable   affine    coordinate 

system   in  the  ^,  ,  y_   -plane   our  mapping   is    then  described 

by   9*1  =  k]?l'    92  ~  k2f  2'    S0   th£it  Fl'   F2  are  of  ths   form 

Fl(fl^    =klfl+    •", 

P2(fl'f2)    =   k2?2   +    •••     . 
We   distinguish   three    cases   according   to    the    signs   of  the 
eigen  values  k,  ,    k0.      In   the    case    (i)    of    two    negative   eigen- 
values,   k„  <  k,    <  0,    every    solution  f  which  comes   sufficient- 
ly  near   to  the    origin  plunges    into   it  with  x  tending    to 
+  ^^  ;    its   asymptotic   behavior   is  described  by 

<Z\(x)-a.cklx,         c72(x)-J.eklx        (x  ->  o©) 


in  the    sense    that 


-k,  x 


-v 


+   0(e^X), 
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where  £    is    a    positive    cc  ,  for   all   solutions, 

e    constant   a    is    specific    ( node ,    see    Pig.    5).      In 
the   case    (ii)    of  one    negative,    one    positive   eigen-value, 
:•:-,    <  0,    kg   >  0,    there    exist   exactly    tv.o    solutions 
(1  +  ,    1_   in   Pi  '.    6)   appro  ag  0   with  x  — -»   +<&}    their 

asymptotic  behavior   is   described   by 

fY ( :•:  )/nh-    e V ,        ^(x)-vO.e ^        (x  -*.   +  <^»  ) 
( saddle ) .      Of    course    this    statement    is  to    be    interpreted 
so    that   equivalent    solutions   are    considered    as    one    and 
the       .      .  addle    is  better  understood   by  observ- 

ing  that    in    :■]  L       ■    se    (ii)    two   other    solutions,    2      and   2    , 
plunge    into    0   in   the   >   Lreetion  of    the  /^-axis,    while   x 
goes    to    -  <yr.  se    (iii)    of    two    positive   eigen-values 

there    is   no    solution  whatsoever   w]   i  ■      approaches 'the    origin 
as    x    tend;:    to    +  <:■& ,     ;xcept    the    trivial   one   ^C.    =  V?    =   0 

ore  as,    of    course,    every  solution  comes    sufficiently 

near   0    tahe  s  lunge    when  x   moves    toward    -  oo), 

V.'ith  .•        in  mind    we    now    advance    our  de- 

cisive 

6 .  ic    different! a  1_  equations    of 

the   s.^ocl:   layer  del  c     ..  ition  from   a    state    Z 

to    a    state    Z-^   >  ^c,    the   initial   state    ZQ    Is   a    node,    the 
end    state    Z-,    a    saddle. 

Cle  behavior    in   the   neighborhood   of   the 

two    singular  Z    ,    Z-,    is   sue]  to   favor   the  existence 

of    a    unique    solution  of    J  Layer   pro    Lem.      All   other 
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combinations   would   make    one    expect   either    no    solution   or    an 

infinity  of   solutions. 

Again,    the   proof   depends,    above    all,    on    the   funda- 

o 
mental   inequalities    (14)   for  M      =    (p,    -   P0)/(v0    ~   vn  )   which, 

7,./'V      being    the   adiabatic  derivative   m   =   -  -M  we   may  write 

in    the   form 

(35)        m27(°>-   7<°>  >0,      M2,!1'  -  7«d  <0. 

They  must    be    combined   with 

(56)  I.    >  Cj  Jv  >  0,       Tp   >   0 

First  investigate  the  neighborhood  of  Z„  and 

°  o 

therefore   write 

;;   (V    -   a)    +  p    =  M2(V    -   VQ )    +    (p    -    pQ). 
For   infinit'  ]    &\    =  p  -  p    .        Sv  =  V   -  V 

$ E   =   T0^/        -  pQ  <Tv,    <fi(V   -  a)2   =   CV0    -  a)^V,    hence 

therefore    in   neglecting    terms    of   higher  order 

E    -  ^H2(Y   -   a)2    -    c    —    To(^-    yQ). 
The    following   approximate    linear   system  for  <f  V,     <3p  results: 

fe(/v)  =  ^0U:E<Tv  +  Ip), 


ere 


^   =   1/tyu*   >   0,       ^'    =    TM/£     >   0 


and  cPv,     J*  I   stand    for 

(o)     *  „(o)   C  c  (0),  (o)   C 

^     =    ^p         SP   +    7V       *V,         >T   =    T^     £P    +   Tv       ^  V. 

In   order    to    avoid    unnecessary  encumbran  ce  ,   we   drop   the    index   0, 

Then    the    equations   determining   the    eigen-values   k  and   corres- 
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j  one 


ponding  eigen-directions    (<£  V,     S  p)   at    Z      Decora 

k.    f  V  =   *(M2  <f  V  +     Sp),        k.     <Tt   =  ^».  /7 

giving   rise    to   the    following   secular  equation  for   k, 


k   -    fM    , 


r 


=  o 


or 


-     / 

kTv  -   f'lv'     kTp"   f  7p 
(37)    Tp.k'    -    [^  (H2Tp    -    Tv)    +    ^'    7p).k   +  f  flvfyv    -   %)    =   ° 

Compute   the  discriminant 

A=  {^»2Tp  -  V  +,f  7p)e  -  4ff  V"87p  -  7y>  - 

{f<:i%   "   V    "f  7p}2   +  Vf  <TP    ?v   "   Tv    7p» 

and  thus  verify  by  (29)  that  /^   is  positive, 

A  -{  yte2Tp  -  Tv,   -   f'   y^*  +  4ff   s  iff   . «,  >  0. 

Consequently  ■  the  eigenvalues    are   real  and   distinct,    at    the 
singular   point  Z      as   well   as   at    Z-,* 
Their   product   equals 

W"%    ~     7v}/Tp, 
and   because    of    (36)    and   the   fundamental   inequalities    (35) 
that    product   is    positive    at   Z    ,    negative    at    Z-,.      Hence    the 
two   eigen   values    are    of    opposite    sign   at   Z-,,    of   equal   sign 
at    Z    .      V/hether    the    latter    sign   is    positive    or  negative    can 
be    decided  by   the   sum   of    the    two    k's   for   which    (37)    gives    the 


Value 


fy(M2Tp   -   Ty)    +    ff'%)    /T3 


All   constituents   of   this   ex  Loi    are    positive,    including 

o    (0)  (0)         viLo)  (o)  (o)  iA  (o) 

Thus   we    arrive    at    the   result    that   the    two  roots    k  at   Z     are 


positive.      Our   theorem  has    now  been  proved,    and    a    third   aspect 
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of    the   fundamental   inequalities    (15)   has   been   disclosed;        lot 
only  do    they  rule    that   by   a    shock   the    flow    changes   from   super- 
sonic   to   subsonic   and    that    glancing    incidence    of  a    shock  wave 
reflected   by  a   wall   is    forbidden,    but    they  also    decide    that 
one   of   the    two    points    Z    ,    Z,    is    a   node,    the    other    a    saddle. 

In  order  actually   to    prove    the   existence    of    a    unique 
solution  of  the    shock   layer    problem,    one    has    to   bridge    the 
gap   between    the   neighborhoods    of    the    two   singular   points    Z    , 
Z-,    by   some    sort    of   topological   argument.      In   that   respect   one 
fact    is   of  decisive    importance: 

■ho  or  em  7 .      _hhj_  d  if  fere  ntial  equat  ions    of   the    shock 
layer  have    no    sir.guhr  po int   besides    Z   ,    Z -,  . 

Indeed    if    there   existed   three    distinct    singular 
points    Z   ,    Z-,,    Zp    they  would    satisfy    the   equations 
7o    +  p0    =  I:I2V1    +   Pl   =  I.I2v|    +  p2   =  K2a, 

Eo    "  l:i2(Vo    "   a)2    =  31    "  ^2(V1    "    o)2    =  E2    "  Im2(V2    "    a)2' 

The  first  set  states  that  Z  ,  Z-,,  ZQ  lie  on  a  straight  line 

2 

of   negative    inclination   -  M    ,    and   we   may  assume    the   arrange- 
ment   Z      <  Z-j    <  Zp.      Then   the   second    set   requires    that    both 
Z    =  Z,    and    Z?    satisfy   the    Hugoniot    equation   H(Z    Z    )    =0, 
in   contradiction    to    a   st  pa  jviously   proved   by  vhioh 

there    cannot   be   more    than    one    point    Z   on   any  ray   from  Z 
satisfying   that    equation.         The  observation   is    urged    upon   us 
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that   whereas    the-   Hugoniot   relation  is    not   transitive,    the 

relation    between    Z    ,    Z-,    defined    by   the    simultaneous    equation: 

JV0    +   P0    =  M2V1    +   P]_,  H(Z1Z0)    =   0 

is.(M      being    a  given    constant). 

The    topological  picture   of  the    trajectories  which 
one    expects    is    indicated   by  Pig.    7.      If    this    is    correct,    the 
re  -ion   covered   by  all   tra  iectories   running    into   Z     for 
x  —a,   -  o^w°uld   be   bounded   by   the    saddle    line    through   Z-, 
consisting    of    the    two   branches    2      and    2_    (Fig.    6).       The    out- 
look is   encouraging. 
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